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HYDRODYNAMIC EMPACT OF A SYSTEM WITH A SINGLE ELASTIC MODE 
I -THEORY AND GENERALIZED SOLUTION WITH AN APPLICATION 

TO AN ELASTIC AIRFRAME ^ 

By WiLBTiR L. Mato 


SUMMARY 

Solutivns of impact of a rigid prismatic float connected by a 
massless spring to a rigid upper mass are presented. The 
solutions are based on hydrodynamic theory which has been 
experimentally confirmed for a rigid structure. 

Eguations are gicenfor defining the spring constant and the 
ratio of the sprung mass to the lower mass so that the two-mass 
system protides representation of the fundamental mode of an 
airplane wing. The forces calculated are more accurate than 
the forces which would be predicted for a rigid airframe since 
the effect of the fundamental mode on the hydrodynamic force 
is taken into account. The response of the two-mass system 
gices the response of the represented mode and, although no 
proriswn is made for taking into account the effect of secondary 
modes on the hydrodynamic force, means are indicated whereby 
the results may be used to approximate the response of modes 
other than the fundamental mode. 

Time histories of the hydrodynamic force and structural 
response are giren for wide ranges of mass distribution and 
ratio of natural period to the period of the impact. By use of 
nondimensional coefficients these results are made applicable 
to different combinations of velocity, weight, angle of dead rise, 
and fiuid density. Although the equations permit solutions for 
different combinations of fiight-path angle and trim, an approxi- 
mation is gicenfor correcting the results for the combination for 
which solutions are giren to other conditions unihin a narrow 
range indicated to be of primary interest to the design engineer. 

In a comparison of the theoretical data with data for a severe 
flight-test landing impact, the effect of the fundamental mode on 
the hydrodynamic force is considered and response data are 
compared with experimental data. Consideration of the funda- 
mental mode alone fails to account for the fact that during the 
impact partial failure of the inboard-engine mounts occurred, 
but use of the theoretical solutions to approximate the effects of 
further wing torsion leads to substantial agreement. 

INTRODUCTION 

In recent years the development of large airplanes has 
caused the elastic behavior of airframe structures during 
landing impact to become important. The work which has 
been done on this problem has been handicapped by lack of 
proper knowledge of the time history of applied ground 
reaction. This situation has been particularly acute for 


seaplanes because of difiSculties in measuring the hydro- 
djTiamic force, the seaway, and the manner of contact with 
the seaway'. 

In order to facilitate the interpretation of Si^t data and 
to lead to the prediction of design loads on a rational basis, 
a theoretical hydrodjTiamic studj" was made and tests of a 
rigid float were conducted at the Langley impact basin. 
Since the results of these tests agree with the theoretical 
results for wide ranges of the pertinent variables in numerous 
force time histories (reference 1), it is assumed that the 
theory may also be used in considering the effect of the 
upper-structure elasticity of a seaplane on the motion and 
force characteristics of the hull proper, which is assumed to 
be rigid. 

The bending of wings during impact, which for modem 
flying boats is tlie primary structural action, is considered in 
the present report by reducing the fundamental mode to an 
equivalent two-mass system. The results are presented in 
a form suited to general application and are compared with 
experimental results for a particular case. The equations 
showing the method of solution are included in appendix A, 
and a sample data sheet is given as table I. 

SYMBOLS _ . 

t„ time required for one-fourtb cycle of natural 
vibration 

tf time between initial contact and maximum hydro- 
dynamic force for rigid structure 
t time elapsed after initial contact 
iWi lower, or hull, mass of two-mass system 
mg upper, or spnmg, mass of two-mass system 
m gross mass (JVIS or mg+niL) 

W gross weight 

g acceleration due to gravity 

K spring constant of spring connecting mg and wil, 

force per unit deflection 

n, acceleration normal to water surface of nodal point 
of elastic system, multiples of acceleration of 
gravity; for two-mass system, acceleration of 
center of gravity 

n» oscillatory acceleration of hull about center of gravity 
of two-mass system or nodal point of represented 
mode, multiples of the acceleration of gravity 


‘ Sup«cs«<les NACA TN 1398, ‘‘SoIatbHU Ibr Hfdrodjmaiiilc Impact Force tad Beaponse ol t Tto-Mies Sj^stem. irlth an Application to an Elastic Airframe'’ by Wilbur L. Mayo, 19{7. 
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resultant velocity at instant of contact with water 
surface 


P 

T 

7o 




C, 

c, 

c, 

y^nax 

/» 


mass density of fluid 

angle of trim; angle of hull keel with respect to plane 
of water surface 

flight-path angle at contact; angle between flight 
path and plane of water surface 
angle of dead rise 


nondimensional time coelRcient 



nondimensional load-factor coefficient 



nondimensional draft coefficient 

f 



draft at instant of maximum acceleraticJn . 
natural bending frequency 


Where units are not given, any consistent system of units 
may be used. 


THEORY 

HTDRODVNAMIC 

The hydrodynamic theoiy used in the present report is the 
same as that developed in references 1 and 2.. A basic 
differential equation which gives the instantaneous force in 
tenns of tlic instantaneous position and motion of the float 
is given in reference 2. This equation is used herein to 
determine the effect of airframe elasticity in. altering the 
motion and force time history (appendbe A). The solution 
is based on the assumption that the float does not change 
trim during impact. In this comiection .the pitchir^ moment 
may be large, but the time of the impact is short enough to 
warrant (at the present stage) neglect of the resulting angular 
velocities and displacements. 

The solution presented herein is for a prismatic float 
with such beam loading that the chines do not immerse^ 
during impact. For waves that give the severe design 
condition of full-length impact, conventional beam loadings 
arc small enough to cause the maximum force to occur at 
drafts sufficiently small to make the effects of finite width 
and chine flare secondary. Reference 1 indicates that for a 
conventional float neglect of the pulled-up bow is justified 
when the trim is 3° or. greater. Although for high-trim 
landings initial contact by the afterbody may substantially 
change the trim before the . main forebody impact, the 
neglect of afterbody loads is justified because, during the 
main impact, the shielding of the afterbody by the forebody 
due to depth of the step and to keel angle is such as to mini- 
mize the importance of afterbody loads. 

STBUCTCRAL 

A simplified representation of primary elasticity of an 
airframe is shown in figure 1. A rigid lower mass is 
considered to be connected by a massless spring to a rigid 
upper mass m^. In determining the fundamental bending of 
airplane wings part of the wing mass must be included in 
mi and part of U)e wing lift should be applied to m^. In 
the present report the gravity force on each mass is assumed 
to be balanced by wing lift. 



Fiouae 1.— simplified representation of primary elastldty of an elrtnuno. 


The problem of determining the properties of the two- 
mass system so ’ that it is representative of the primary 
elastic action of the airplane is rather simple if it is assumed 
- that during the impact the structure deflects with the shape 
of its fundamental mode of vibration. The requirements arc: 

(1) The total mass of the simplified system must equal 
the total mass of the airplane in order that the proper nodal 
or center-of-gravit}' accelerations can bo obtained. 

(2) The energy of vibration for the same amplitude of 
the hull and lower mass (relative to the nodal poini) must 
be the same for the two-mass system as for the considered, 
mode of the airplane structure. 

(3) The natural frequency of the two-mass system must 
be the same as the frequency of the considered inode of 
airplane vibration. 

Equations which permit determination of the masses and 
sprir^ constant of the simplified system so that it meets 
these requirements are given in appendix B. These equa- 
tions and the foregoing requirements are applicable for both 
landplanes and seaplanes. 

In the present report the represented structural mode is 
considered to bo devoid of vibration prior to the instant of 
impact. Thus, the computations may repnaent cither a 
first impact or a subsequent impact resulting from a bounce 
• sufficiently high to cause aerodymamic and structural 
damping to stop the vibration during the time the seaplane 
is in the air. This report does not give a representation of 
successive impacts, such as might occur in seaways, which 
lead to accumulative or resonant effects. Available flight 
data indicate that a single heavy impact, such as that 
considered herein, is the primary' cause of structural failul’e.s. 

The response of the two-mass- system is obtained in <-on- 
nection with the calculation of the time liistory of the hydni- 
dymamic force, and from this result the complete rosoonsc 
of the represente-d mode can be obtainctl by' the simple 
procedure given in appendix B and demonstrated in tin* 
section entitled “Comparison with Experiment.’’ The 
response of other modes to the force computed on the basis 
of the fundamental mode can be separately determined 
and superposed (reference 3). In order to minimize the 
complexity of the solution, how'ever, the present investi- 
gation does not provide for taking into account the effect 
of the other modes on the hydrodynamic force. Although 
the other modes may have a substantial effect <jn the local 
loads in the structure, the effect of these modes on the 
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hydrodjuamic force is considered to be secondary as com- 
pared with the effect of the fundamental mode. 

If a large number of solutions for the two-mass sj’stem 
have been made in order to determine the effect of the funda- 
mental mode of different ^\■ings on the hydxod\Tiamic force, 
the response of modes other than the fundamental can be 
appro.ximated from the use of a solution for a mass ratio 
and ratio of the time period of the impact force to the natiual 
period of ihe two-mass system representative of the con- 
sidered mode. An e.xample of such use to approximate the 
effects of wing torsion is given in the present investigation 
in a comparison of results of computations with e.xperi- 
mental results. 

If the response of more than one mode is considered, the 
structural and aerodynamic damping, which are not con- 
sidered herein, are important factors in determining the ex- 
tent to which the maximum response of the different modes 
should be superposed without regard to phase relationship. 
It is expected that the effect of the damping will be most im- 
portant for the higher modes and that a result leading to con- 
servative design will be obtained if damping is not considered 
and the maximunLs of the first two or three modes are super- 
posed without regard to phase relationsliip. 

RESULTS 

Solutions of the equations in appendix A were made for 
wide ranges of the pertinent variables. Time histories of 
the calculated nodal acceleration, or hydrodynamic force in 


terms of the weight, are given in figures 2, 3, 4, and 5 for 
ratios of the sprung mass to the lower mass equal to 0.25, 
0.60, 1.00, and 1.36, respectively. Each figure is three- 
dimensional; the third dimension is tjtt which is a ratio of 
the period of natural vibration to the speed of the impact. 
In representing the period of natural vibration, f, is taken as 
the time required for one-fourth of a cycle. The speed of 
the impact is represented by making tt equal to the time be- 
tween initial contact and ma.ximum acceleration for a' rigid 
structure. If the time to reach ma.ximum force for the elastic 
structure should be used in defining tt, discontinuities in the 
time to reach ma.ximum force would cause discontinuities in 
the time-ratio scales of the plots. (See figs. 2 to 6.) 

An expression for f, may be obtained from the relation 

and equation (A4) in appendix A. The equation 
for tt is as follows: 




( 1 ) 


where 

time coefficient at instant of maximum acceleration 
for rigid body (0.678 for ^=22.5°, 70= 14°, and 
t=3°) 

The e.xpreseions for and tt may be used to determine that 


V I'll 

/ mam A 

tt 2C. 





^ 2 ) 



352 


HEPOHT 1074 — NATIONAI- ADVISORY COMMITTEE FOR AERONAUTICS 



The oscillatoi'y acceleration is given in figures 6 to 9. 
These figures are the same as tlie figures giving the nodal 
acceleration (figs. 2 to 5), except that the acceleration plotted 
is the diffei'cnce betiveen the hull acceleration_and the nodal 
acceleration.. The time Iiistoiy of the.huU acceleration can 
be obtained by summing the two plots for a particular mass 
ratio. On the basis of the assumption that during impact 
the airframe structure deflects in a particular mode, a time 
liistory of the acceleration of any point in the structure is 
obtained from the results for the equivalent two-mass system 
by the following procedure: 

(1) From the deflection curve of the represented mode, 
obtain the ratio of the deflection of the point of interest to the 
deflection of the hull. Both deflections are taken relative 
to the nodal point. 

(2) Multiply this ratio by the oscillatory acceleration 
given either by figure 6, 7, 8, or 9 or by interpolation between 
these figures for the mass ratio of the equivalent two-mass 
system. 

(3) Add result to the nodal acceleration given by figures 
2 to 5. 


Time histories of the acceleration given in figures 2 to 9 
are on a nondimcnsional basis. The nondimensional cQr 
efficients, wliich contain velocity, weight, fluid density, and 
acceleration of gravity, were used in reference 1 in . a com- 
parison of theoretical data with impact data for a float 
having an angle of dead rise of 22^°. 

Application to other angles of dead rise.— The function of 
the angle of dead rise can also be included in. the nondimen- 
sional coefficients, but in the present investigation thb 
function is isolated and treated as a factor for correcting tlie 
results presented for angle of dead rise of 22J° to other dead- 
rise angles. The pertinent relationships between results for 
dillereut angles of dead rise may bo expressed as follows: 


r 1 1*'® 

(3) 

n,ccm)f{AW 

(4) 


w'here 

/(jS) function representing variation of virtual mass for 
two-dimensional flow with angle of dead nse 
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f(A) aspect-ratio factor including effect of angle of dead 
rise on aspect ratio; ratio of Tirtual mass for 
three-dimensional flow to virtual mass for two- 
dimensional flow 

acceleration at any point, either oscillatory, nodal, 
or total 

The shape of the force and acceleration curves for a given 
value of tjtt is independent of angle of dead rise, but the 
effect of angle of dead rise on fj, as given by relation (3), 
does enter into the determination of the value of tJtt for a 
particular solution. After the value of tjft for a particular 
solution has been determined, the acceleration and time 
values for an angle of dead rise of 22J° are proportioned by 
means of relations (3) and (4) to the corresponding values 
for the angle of dead rise used in determining f,/f« in order to 
obtain the proper acceleration history. , 

Although adequate impact data have not been available 
for checking the theoretical equations for angles of dead rise 
other than 22\°, the theory is equally applicable to planing 
floats. Study of planing data has shown that the functions of 


angles of dead rise used in equations herein are approximately 
correct for angles of dead rise ranging from 15° to 30°. The 
functions are: 



( 5 ) 


( 6 ) 


Until improved functions of angles of dead rise are obtained, 
functions (5) and (6) should be substituted in equations 
(1) and (2) and in relation (4) to correct for angles of dead 
rise within the range from 15° to 30°. Eough approximation 
can be obtained by use of functions (5) 'and (6) for angles of 
dead rise greater than 30° but not for angles of dead rise much 
less than 15°. In reference 2 there is a discussion of the 
inadequacy of f(A) for aspect ratios which normally occur 
for small dead-rise angles. For very large angles of dead rise 
and moderate velocity the static forces, which are not 
considered in the present investigation, become of greater 
importance. 


2724SS— 54 ^24 
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Approximate correction to other flight paths and trims. — 
The flight-path and trim conditions of primary interest to the 
designer tend to be independent of the variables in seaplane 
design and constitute a limited range within which approx- 
imate correction of results to different. angles of flight path 
and trims can be made without necessity for repeated time- 
history solutions. For any particular combination of 
horizontal speed, rate of descent, and tiim, the most severe 
impact load for most of the structure occurs when the sea- 
way is such that the keel contacts a wave slope approximately 
parallel to it. The effective angle of flight path and trim for 
such an impact are defined relative to the inclined wave slope j 
therefore, the trim which gives maximum force is zero. The 
largest flight-path angle, relative to . the keel and to the 
critical wave slope, is also associated with the most severe 
force. The value of the largest flight-path angle is not so 
definite as the critical trim. but tends to be independent of 
variations in size and wing loading. The velocity of the wave 
should be considered in determining the contact speed and 
flight-path angle. 

The equations and method of solution given in appendix 
A permit solution for different flight-path angles and trims; 
however, approximate correction of the results in figures 2 to 
9 to other positive contact angles can be made by assuming 


tha,t the proportionate effect of the sU'uctural elasticity on 
the hydrodynamic force is solely dependent on the ratio 
tJU. Curves given in reference 1 show values of (7. f or 
different flight-path angles and trims, which may bo sub- 
stituted in equation (2) to obtain the value of tJU for differ- 
ent contact angles. In making the approximate correction, 
the solution presented herein for the obtained value of the 
ratio tjU should be used to approximate the shape of the 
curve giving the desired time historju The load or accelera- 
tion scale should be corrected to the different contact angles 
by multiplying the present result by the ratio, determined 
from curves given in both references 1 and 2, of the load- 
factor coefficient for the different contact condition to the 
load-factor coefficient for the conditions of trim and 14° 
flight-path angle considered herein. Correction of the time 
scale involves a similar procedure in which the time coef- 
ficient is used rather than the load-factor coefficient. 

Since the force curves for a rigid body are approximately 
the same shape for different angles of flight jiath and trims 
(reference 1), the approximate method of corrcctuig to dif- 
ferent angles of flight path and trims would be almost correct 
if the sti'uctural elasticity did not affect tlie hydrodynamic- 
force curve. The . percentage change in the force on the 
float due to elasticity is a function of the percentage change 
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in draft caused by the elastic compression for a given center- 
of-gravity position. An indication of the validity of tbe 
approximate correction is obtained by studying the extent 
to which the ratio of the spring deflection to the draft is 
constant for impacts of the same values of tjti at different 
angles of flight path and trims. From the expressions for 
tjtt, Cl, Ct, and Ct with the spring deflection assumed to be 
proportional to the hydrodynamic force, this ratio may be 
represented by the expression CiOtlCt, in which all values 
are for the instant of maximum acceleration. The variation 
of this expression with flight path for an impact of a rigid 
float at 3° and 12° trim is given in figure 10. Values of Ct 
and C\ used in obtaining this figure are given in reference 1 ; 
values of Ct were obtained in conjunction with the data of 
reference 1 but have not been published. 

In the present report the numerical values of CiC*(Ct 
have no significance and they are of interest only because of 
the extent to which they are constant. Figure 10 indicates 
that for large flight-path angles and small trims the ratio is 
approximately constant. The deviation from a constant 
value of this ratio is due to planing forces which exist in an 
oblique impact and become more important for low flight- 
path angles and high trims. The conditions of laige flight- 
path angle and small trim previously adjudged to be of 
primary interest to the designer constitute the ranges in 
which the deflection ratio is fairly constant and the appro.xi- 
mate method of correction should give a fair degree of 
accuracy. The present solutions are considered to be for 
conditions suited to correction of the results to other con- 
ditions of greatest practical interest; they represent a moder- 
ately severe combination of flight-path angle, w'ave slope, 
and trim chosen to facilitate correlation of the theory with 
an impact which resulted in substantial damage to a well- 
instrumented flying boat during flight tests. 

The equations presented herein are not valid for zero 
trim; an assumption that the float is prismatic gives solu- 
tions of infinite w'etted length and infinite force for zero 
trim. Solution for 3° trim and a prismatic float is much 
simpler than a correct solution for 0° trim because necessity 
for consideration of bow' shape is eliminated. The solution 
for 3° trim may be taken as an approximation of the critical 
design load or, as illustrated in the following section in a 
comparison of theory wdth experiment, an empirical factor, 
which includes bow effects, may be used to convert values 
of acceleration and time for 3° trim to values for 0° trim. 

COMPARISON WITH EXPERIMENT 

Impact basin data have not been obtained for suitable 
models, and most fli^t landing data have been inadequate 
for the present study. The only data whicli appeared suit- 
able for this comparison are those which W'ere obtained witli 
a four-engine fljing boat, the data for which have not been 
published. Data were obtained for a large number of test 
landings, but only one of the impacts is very weU suited to 
tlie present analysis. This impact gave loads sufficiently 
big b to cause large effects of elasticit}* of the W'ings. The 
impact occurred against the flank of a sizable wave (4 ft) 



Flrghi-poth ang/e, ri, deg 

Fiocxe 10.— Vaiifttkia of lbrce.I)enetrstlan ratio at maitmum fiirce wltli fllgbt.patli angis 

and tiim. 

and thus facilitated the use of results based on a planar water 
surface. The present comparison of theory with experiment 
will be restricted to thin impact since other impacts involved 
more complex contact with seaway and gave less force. 

The horizontal speed, rate of descent, and trim were re- 
corded. A large number of pressure instruments distributed 
in the hull permitted determination of the water surface 
relative to the hull. Data recorded by these instruments 
indicated that the wave slope in contact with the hull w'as 
approximately planar, that the trim relative to the wave 
slope was 0°, that the resultant velocity, considering the speed 
of the wave, was 85 feet per second, and that the flight-path 
angle relative to the wave slope was 14°. 

Structural data available for the test flying boat are not 
adequate for the present anal 3 'sis. Since this flj'ing boat has 
the same number of engines and appro.ximately the same 
gross weight and horsepower as a landplane for which a 
large amoimt of structural data is available, assumption is 
made that the flj ing-boat wing has the same mode shape and 
mass distribution as the \ving of this landplane. Use in 
equation (B6) in appendix B of data for the landplane given 
in reference 3 leads to the following mass ratio of the two- 
mass system representing the fundamental mode: 

^^ 0.25 

niL 

Based on studj’ of the acederations at several points in the 
test flj'ing boat dining periods of rdativdy free vibration in 
which the fundamental wing bending mode appeared to be 
predominant, a natural frequency of 3.6 cydes per second 
was selected for use in the present example. The funda- 
mental mode frequency' of the wing of the landplane is 3.4 
cydes per second. (See reference 3.) 

Since the mass ratio is equal to 0.25, figures 2 and 6 are 
used to approximate the action of the fundamental mode. 
Further, substitution of values for conditions for this impact 
in equation (2) results in use of the specific time history 
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given for ^=1.2. Substitution of tlio contact conditions in 

the load and time coefficients fixes the. load and time scales 
in an absolute sense. Correction from 3° trim to 0° trim may 
be made by assuming that the shape of the time history is 
approximatelj’' the same for both conditions. Empirical 
coiTection of the curves from 3° trim to 0° trim may be seen 
from reference 2 to require a 10-percent reduction in the 
acceleration. An analysis of data obtained at the Langley 
impact basin for impact at 0° triiri indicates that correction 
of the time scale from 3° trim to 0° trim requires a 10-percent 
reduction in the time values. 

Hydrodynamic force. — ^The nodal-point acceleration lit, 
obtained by the procedure discussed herein, represents the 
hydrodjTiamic force, in multiples of the Aveight, applied to 
the flying boat. Since tlie experimental data do not provide 
measurement of the hydrodynamic force as such, direct co.m- 
parison of the tlieoretical force-curve results witli experi- 
mental results is not permitted. Instead, a comparison of 
the theoretical response of tlie structure with the experi- 
mental response is necessarj*, and, if the agreement is ade- 
quate, it may be concluded that both the hydrodynamic and 
the structural actions are adequately represented. 

Before a study of the response of the structure is made, the 
theoretical effect of the response on the hydrodynamic force 
should be observed. This observation is made by comparing 
the force curve obtained for a mass ratio of 0.25 Avith the 
force curve for a rigid structure. Both curves arc included 
in figure 11. Tlie curve, for the case of a rigid structure 

( !H5=o jAA’as obtained from reference 1 for 3° trim and cor- 

\mL J 

reeled to the conditions of the present example as preA'iously 
indicated. Comparison of the curves for mass ratios of 0 
and 0.25 shoAVS that in the present example the theoretical 
effect of the structural elasticity on the maximum hydro- 
dynamic force is to reduce it 15 percent. For the hypo- 
thetical condition of a concentrated wing mass located at a 
point in each semispan of a massless wing structure, the 
conditions of the present example Avould give theoretical 
reduction in the ma.ximum hydrodynamic foi'ce due to 
structural elasticity of 44 percent. This result is indicated 
by comparison of the maximum of the cun-e in figure 1 1 for a 
mass ratio of 1.36 AAnth the maximum of the cuta'c for mass 
ratio of 0. 

The curves in figure 11 sIioav only reduction of the hydro- 
djmamic force because of structural elasticity. It should 
not be concluded, however, that the effect ia alAvays in this 
direction. Figures 2 to 6 shoAv that in some cases the hydro- ' 
dynamic force is increased; the maximum mcrease which Avas 
calculated Avas of the order of 12 percent. 

Hull acceleration. — By combining the results for the two- 
mass case with approximation of the pitching action on the 
basis of a rigid structure, the following equation may be 
obtained for the hull acceleration at different longitudinal 
stations: 

( 7 ) 



0 .04 .OB .iS .tG .80 

Time, sec 


Figu&s 11. — Variation of hydrodynamic force with time. r»0®; T'i*88 fL-et per 

second; Tr-40,000 pounds; cycles per second. 

where 

n* hull acceleration 

/=l+g 

's distance from station to center of gravity 

k radius of gyration (12 ft in present exampli') 

I distance from resultant hydrodjmamic force to 

center of gravity. 

The forebody length of the flyhig boat is 31.75 feet. Be- 
cause of hoAV effects, a length of 25 feet is assumed to have a 
rectangular loading for this zero-trim impact and the result- 
ant force is located 12.5 feet foiward of the step which leads 
to a value of I equal to 8 feet. Use of the foregoing procedure 
to calculate time histories of the acceleration for two stations, 
in the hull at Avhich accelerometers Avere located gives tlnj 
curves in figure 12. The maximum accelerations recorded at 
those stations are also shoAvn; agreement with the computed 
maximum acceleration is good. The full experimental time 
history is not included because tlie film speed Avas not great 
enough to permit accurate determination of the shape of the 
time history. This factor, together with some uncertainty. _ 
in defining the exact instant of contact, prevents exact check 
of the time to reach maximum acceleration; thus, the experi- 
mental points in figure 12 are located at the theoretical 
time of maximum acceleration. 

Although the agreement of calculation Avith experiment in 
figure 12 is good, this agreement can be interpreted as con- 
firmation of the elastic action of the structure, only to the. 
extent that disagreement of hull accelerations computed on 
the basis of a rigid structure can be shown. Values of 
maximum accelerations computed for a rigid structure are 
included in figure 12. Tlie disagreement with experiment 


n^=nj+n, 
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Figuec 12.— Varlatfon of bow md step accelcratloD with Hine. 


is greater than for computations in which the elasticity of the 
structure is considered. The difference is small, however, 
because the effect of elasticity in increasing the hull accelera- 
tion for a given force is largely offset by the effect of elasticity 
in reducing the hydrod\"namic force. Since the probable 
accuracy of available data must be considered, experimental 
proof of the theoretical effect of the elasticity on the hull 
acceleration cannot be claimed. Support for a conclusion in 
this respect, however, can be obtained by conjparing the 
theoretical response of the wing with the experimental 
response. .If the wing responds as assumed, the basic equa- 
tions require that the hydrodjuamic force and hull accelera- 
tions be as calculated. 

Elastic axis. — Use of the nodal and oscillatory accelera- 
tions of the representative two-mass system to predict 
accelerations along the elastic axis of the wing requires 
consideration of the fact that the elastic axis of the wing is 
mu at the center of gravity of the flying boat. An approxi- 
mate correction may be obtained by midtiphdng the results 
for the two-mass case by the factor/ which is used in equation 
f7). This correction is not entirely consistent with that 
given by equation (7), but each approximation is considered 
more accurate for its particular case. An improvement to 
the present correction which would change the results 2 to 3 
pen.‘ent might be made, but the complication is not considered 
to be warranted. 

In the present example the value of/ is 0.S6. Application 
of this factor, of equation (B7) of appendix B, and of per- 
tinent stnictural data given in reference 3 gives the curves 
in figure 13 in representation of the acceleration time histories 
of the elastic axis of the wing for the hull, inboard-engine, 
nodal, outboard-engine, and tip stations. 

Wing torsion. — The torsion of the wing during impact may 
have substantial effect on the acceleration of engine and 
nacelle masses forward of the wing. Use of the procedure 
which gave the acceleration time histories in figure 13 to 
calculate the acceleration at the engine gives a maximum 
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acceleration of Zg at the inboard engines and a maximum 
acceleration of 5.6p at the outboard engines. The relative 
magnitude of these values is in strong disagreement with 
the fact that during this impact partial failure of the inboard- 
engine mounts occurred, but the outboard-engine mounts 
were not damaged. 

Since the structural data used in the preceding computa- 
tion are for the actual fundamental mode, their use in impact 
calculations involves assumption that the coupling between 
the torsion and bending corresponds to the coupling which 
exists in natural vibration. Actually, the torsional defiec- 
tions in impact are determined not merely by the bending 
deflections but also by the large nodal acceleration, which 
does not exist in natural vibration. 

In the present example the results which have been calcu- 
lated for impact of a two-mass elastic sj’stem will be used to_ 
predict the response of the engines. The procedure for 
doing this is to select the proper solution and then to adjust 
the acceleration and time scales of the two-mass solution so, 
that the maximum acceleration of the nodal point corre- 
sponds to the maximum acceleration of the elastic axis at the 
engine station. The acceleration of the upper mass of the 
two-mass system then represents the response of the engines; 
however, because of the eccentricity of the impact, an in- 
j crement must be added. 

Data obtained from Guggenheim Aeronautical Labora- 
toiy, California Institute of Technology include the torsional 
;■ deflection of a station inboard of the outboard engine for a 
given moment applied at the wing tip. For a flying boat 
the absence of cut-out for the landing gear tends to give a 
, stiffer wing; therefore, in the . present example, deflection 
measured on the \\ing of the landplane slightly inboard of 
the outboard engine is considered applicable to the outboard- 
j engine station. Relative deflections between the inboanl 
and outboard engines are estimated as follows: 

(1) Torsional deflection at inner engine equal to 1 unit 
due to each engine, or 2 imits total 

(2) Torsional deflection between inboard and outboard 
j engines equal to 2 units because of greater distance of flexure, 

! boosted to 3 units because of increased fie.xibility of structure. 

On the basis of data from California Institute of Tcch- 
nologj' the average of the static moments of the inboard and 
outboard engines is taken as 22,000 foot-pounds. After 
determination of the torsional deflection at the outboard 
engine for this average moment applied at the wing tip, 
multiplication by the ratio 1/2 gives a value for the static 
deflection of the inboard engine and multiplication by the 
ratio 5/4 gives a value for the static deflection of the out- 
board engine. 

An approximation to the calculation of the response of the 
engines to the total acceleration of the elastic a:ds is to 
neglect djuiamic interaction between the engines and treat 
each engine as a single-mass oscillator having a natural fre- 
quency determined by its static deflection. Such a proce- 
dure gives values of 7.6 cycles per second for the inboard 
engine and 5.1 cycles per second for the outboard engine. 
Use of these frequencies to determine and division of by 
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Tims, sec. 

Figxtbs is —V ariation of accelerations at stations along t&e wing elastic ails wltli time. 


ti (t{ equal in this case to the time to reach maximum acceler- 
ation of the curves in figure 13 for the station in question) 
leads to values of tJU for the inboard and outboard engines. 
The next step is to select a mass ratio which for these values 
of fjt{ has a shape of the nodal acceleration time liistory, 
which approxima tes the shape of the acceleration time history 
for the elastic axis at the station in question. In the present 
example the mass ratio of 0.25 is used. Scale factors for both 
the load and time scales axe deterihined so that the maximum 
nodal acceleration for the two-mass solution vdH agree with 
the maximum acceleration and time to reach maximum 
acceleration of the elastic axis at the station in question. 
After these factors are applied to both the nodal and oscilla- 
tory curves for the selected mass ratio and time-period ratio, 
use of the results and equations (B7), (B5), and (B2) in 
appendix B to calculate the acceleration of the sprung mass 
of the two-mass system gives accelerations of the engines. 
Appro.ximation and superposition of the pitching action on 


the basis of a rigid structure leads to the solid-line curves 
given in figure 14 for the accelerations of the engines. 

Also included in figure 14 is the design ultimate acceleration 
for the engine mounts. Comparison of the calculated engine 
accelerations vnth this value shows agreement _of the calcu- 
lation with the fact that partial failure of the inboard-engine 
mounts occurred but the outboard-engine mounts were not 
damaged. 

During the impact an accelerometer was located at the 
outboard-engine station intermediate between the clastic 
axis and the engine. A calculated time history of the acceler- 
ation at the accelerometer location is given in figure 14; this 
time history is based on linear interpolation between the 
computed accelerations at the engine and the clastic a.vis at 
this station in accordance mth the proiwrtionate distances 
involved. The %ure also includes the recorded maximum 
acceleration at tliis point and shows good agreement of the 
computed acceleration therewith. 
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Acceleration time histories for the engmes, computed on 
tlie basis of a rigid structure, are included in figure 14. The 
maximum accelerations computed on this basis do not agree 
with the structural failures which occurred. Furthermore, 
the fact that the curves computed on the basis of a rigid 
structure reach a maximum at the same instant of time is in 
strong disagreement with experiment, which in this respect 
is in approximate agreement with the computations for ah 
elastic structure. 

Difference between the ratio for the inboard-engine 
and outboard-engine stations is primarily responsible for 
difference in the calculated response of the en^es. Most 
of the difference in this ratio for the two stations is not due 
to difference in the sprung-engine frequency but is due to 
the greatlj' different time to reach marimum acceleration tt 
of the elastic a-xis. Agreement with experiment of the 
response calculated at these stations provides indirect con- 
firmation of the acceleration time histories predicted for the 
elastic axis at these locations by the normal-mode method. 

For the impact, e.xperimental data are not available for 
checking the tip acceleration, which is predicted on the 
basis of the normal-mode method, but the initial downward 
acceleration and the 12p maximum acceleration shown in 


figure 13 for this station are in general agreement with 
results recorded in severe impacts of other airplanes. Agre^" 
ment of the computed hull acceleration with experiment has 
already been shown; in an indirect manner all the curves in 
figiu-e 13 exhibit satisfactory agreement with available ex- 
perimental data. Although the response of the engines is 
different from the response assumed in calculating these 
curves, it appears that in practical use the two-mass solu- 
tions given herein can be interpreted on the basis of the 
normal-mode method to obtain both the response of the 
elastic axis and the hydrodynamic force. Further, it appears 
that in practical problems modification of these two results 
is not required when accelerations, of the engine different 
from the accelerations predicted by the normal-mode method 
are determined; the merit of this statement should, be 
independent of whether such modification is made by the 
method used herein or by another method. 

In the foregoing comparison agreement of calculation with 
experiment is obtained without consideration of the response 
of modes higher than the fundamental mode. If, when 
more data are available, it is shown that the response of the 
higher modes can be determined by treatment parallel to 
that given the fundamental mode, the two-mass solution 
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given herein can be used to predict their response by selecting 
a solution, for a mass ratio and natural frequency representa- 
tive of the higher mode, which has a nodal acceleration curve 
of approximately the same shape as the hydrodynamic-force 
curve. 

CONCLUSIONS 

Theoretical solution of hydrodynamic impact of a hull 
mass connected by a spring to an upper mass, the results of 
calculations for wide ranges of mass ratio and natural fre- 
quency, and the use of tliese results in a comparison of 
theoretical data with test data for a flight-tost landing 
impact indicated that: * 

1. In fl>ing-hoat impact the effect of the structural 
response on the hydrodynamic force might ho substantial, 
the shape of the force time liistory might bo considerabli' 
('hanged, and the maximum hydrodynamic force might be 
either reduced or increased. 

2. The greatest reduction in hydrodynamic force occurred 
for the condition of lai-ge mass ratio and low value of spring 
constant. 


3. The normal-mode method was a practical means for 
determining the equivalent twomass system which repre- 
sented the major elastic action of the airframe, for pi'cdicting 
the effect of this action on the hydrodynamic force, and for 
approximating accelerations along the elastic axis of the 
wing. 

4. The acceleration of engines contained in nacelles for- 
ward of the wing could not be computed on the basis of 
coupling between torsion and bending as in the fundamental 
mode, but simple treatment of their response to the com- 
bined translation and oscillation of the elastic axis gave 
agreements with experimental accelerations and gave an 
explanation of partial failure of the inhoard-engine mounts 
during impact. 


LangleyMemorial Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., March 17, 1947. 



HYDRODYNAMIC IMPACT OF A SYSTEM WITH A SINGLE ELASTIC MODE. I — THEORY AND GENERALIZED SOLUTION 363 


APPENDIX A 

MATHEMATICAL EQUATIONS AND METHOD OF SOLUTION 


EQUATIONS OF MOTION 

The following equations of motion, which were derived 
from equadon (30) in reference 2, are for fixed-trim impacts 
of a rigid prismatic float connected by a massless spring to 
a rigid upper mass: 

Acceleration lower mass normal to water surface in feet 
per second per second 

-y,=|sA]/i^y£-t-A:iCosr)*-{-^ [_^+ 

(ms+mi,)[yL—yi^t)+2KiAcosr yi!‘dt+ 

(Al) 

Acceleration of sprung mass normal to water surface in 
fnet per second per second 

2/s=— [(Ay£,®+OTi)yi,+3Ayi-(i/-£,-f-A'tCosr)T (A2) 

Tfl g 

Acceleration of nodal point normal to water surface in feet 
per second per second 


y 




(A3) 


The spring constant in pounds per foot of deflection can be 
expressed by the equation ‘ 


where 

Ki=Xi^smT—yL^ 


rrii^+ms 


(A4) 


cos T 
2 


i/i, velocity of lower mass normal to water smface, ft/sec 
draft normal to water surface, ft 
li angle of dead rise, radians 

T angle of trim, deg 

initial velocity normal to water surface, ft/sec 
initial velocity parallel to water surface', ft/sec 
wis sprung mass, slugs 
niL lower mass, slugs 
/, natural bending frequency, cycles/sec 
p mass density of water, slugs/cu ft 


COMPUTING DIRECTIONS 

A sample data sheet is given as table I. In this table t^ 
numbers in circles refer to rows; the circled numbers under 
the row headings refer to computed values to be used for 
the computations. In the first column the time is equal to ’ 
zero, in the second column the time is Ai, aud in successive 
columns the time is 2 Af, 3 Af, and so forth. Each row is~ 
computed in sequence for any given column before any row 
is computed for the next column, except for the first column, 

TABLE I.— DATA SHEET— GENERAL TERMS 


: Row 


@ 


@ 




Row heading 


' ' ' > L 

Column j 

^ i ^ 

3 



Time, sec 

assumed 

ft 

f = UL* 

cos r)*@X3d 

®X^ 

(mi.+ ms) (©-liijX®) 

^/Ti-4 cos rX@ 

(® +®,) yJ Vi^dt 


SAA'i* cos* rX@ 
®+®+@-(-@ 
(ix®)+® 


(AX®)-t-n»i 

ft/sec* 

C@+@,)f 

@+@, 

@+f 


HIT 

@+^> jIl assumed for next interval 
32.2 


—Ss, g units 


g units 

(mt.X®) + (maX®) 

Tils 


E 


ff units 





364 


REPORT 1074 — NATIONAL ADVISORY COMMITTliB FOB AERONAUTICS 


In the fii’st column the assumed value of is the given value 
of All other values are zero except the values for I'ows 
and @ which in this case are also yt^. Lower-case p 
as a subscript on a row number refers to the computed value 
in the indicated row of the preceding column. 

Eacli row heading indicates the general operation to be 
performed. Wlien the data sheet is set up, numerical values 
sliould be substituted in the row headings for tlie expressions 
that are constant for a specific case. The constants for tlie 
sample computation, given in table II, are as follo\\'s: 


K ^- 107, 664 

Ki — 4.24317 

.4 13a 919- 

ff, deg > - 23.-6 

j/ij,, ft/aec — 20. 6673 

82. 1685 

ma 715.217 

mt 625.776 

fa, cycles/sec - 3 

p, siugs/cu ft - L 938 


For most solutions a time increment At of 0.005 second is 
satisfactory. (The value of At may be varied by considering 


the tune for a given mass to reach a maximum acceleration.) 

■ The number of significant figures to be used should bo 
chosen on the basis of tlie computing equipment available 
and the accurac}' desired. Comparison of a solution com- 
puted with four significant figures with results which had 
been obtained with six significant figures gave a difTereiiee of 
about 1 percent in the maximum acceleration. 

As tlie computations proceed, the lower-mass acceleration 
(row @), tlie sprung-mass acceleration (row @), and tlie 
nodal-point acceleration (row @), all in y units, should be 
plotted against time, (row (i)) in seconds.* Each time a 
column is completed the new points should be added to the 
plot. This plot is tlie only brief method of checking on tlie 
accuracj' of the computation until sufficient solutions are 
obtained to permit cross-plotting. If the points do not lie 
on a smootli curve, an error has been made in tlie computa- 
tions. If an error is made, eveiy value computed thereafter 
contains the eiTor and therefore great care must be taken. 

The number of columns required for tlie computations 
depends on tlie number of columns required for rows @ and 
@ to pass their respective maximums. At least three 
columns should be computed beyond the column in which 
row @ reaches its maximum value. 


TABLE II.— SAMPLE COMPUT.ATION 


[Conatanta are from appendix A] 


Row. 

» 

Row heading 

Column 


1 

2 

3 . 

4 , 

■ 5 

6 

7 

8 

0 

0 

® 

Time At, sec (0.005)* 

— 

(@„xo.ao5)+0, - 

^ + 1 23736“) *© X 40 1 .757 

0X33.4798 

K4O.99[0- (20.6673X©)] 

[(®+®i)0.0025]-t-®„ 

I134.92X® — * 

[(®+®p)0.0025]+@, 

L(@+@.)0.0025H-®, 

7213.62X® - 

®+®+@+® . — 

n50.533X@) + ® 

(133.919X®) +626.776 



(@+@p)0.0025. 

®+®, 

20.6^3+® 

®+f. 

®+f 

r-0. 001398(®-®)]/32.2. 

@/32.2 

(525.776X®) + (715.217X®) 
1240.99 

.,0 

20. 6673 
0 
0 
0 

. 0 
0 
0 
0 
0 
0 

.. 0 

0 

.^0. 

■;o 

b 
- 0 
.,D 
„0 
0 

20. 6673 
20. 6673 

2a 6673 

-JJ 

.,0 

- ^ -0 

a 006 

20. 6673 
0. 103337 
0.010679 
0. 001104 
0. 000114 

2661. 07 
0. 003817 

0 

0. 000003 
0. 003405 
a 000027 
0. 

0 

a 007222 

2662. 16 , 
525. 924 

-5. 06187 
-a 012655 
-0. 012655 
20. 6546 

20. 6483 

20. 6420 

0. 000047 
-0. 157201 

-0. 066575 

0.010 
20. 6420 
0.206579 
0. 042675 
0.008816 
0. 001821 
1061Z 5 
0. 060967 
-0.116663 
0. 000028 
0. 031778 
■ 0. 000160 
0 
0 

-0. 023903 
10608 9 
526. 957. 
-20. 1324 
-0. 062986 
-0. 075641 
20. 5917 

20. 5602 

a). 5287 

-0.000156 
-0. 626230 

-o! 264983 

0.016 
20.5287 
0. 309380 
0. 095716 
01 029613 
0. OO0IB2 
23586.5 
0. 306742 
-0. 781824 
0. 000124 
0. 140730 
0. OOOM6 
0. 000002 
0. 014427 
-0. 319925 
23638 ^ 
629. 742 
—44. 4335 
-0. 161415 
-0.237Q56 
20. 4302 

20. 3455 

20. 268^ 
-0. oo:^ 

-1. 37092 
-0. 5858k 

.. 

0. 020 
20.2688 
0. 411128 
0. 169026 
0. 069491 
0.028670 
40782. 0 
0. 956518 
-2. 75252 
0. 000373 
0. 422190 
0. 001168 
0. 000006 
0. 043282 
-0. 133053 
40581. 7 
535. 082 
-76. 8420 
-0. 300690 
-a 537746 
20. 1296 

19. 9793 

■ 19. 8290 

-0. 008696 
—2. 35534 

-1. 00291 

0. 025 
19. 8290 
0. 611025 
0. 2611147 
0. 133453 
0. 008198 
60767. 4 
2. 28826 
-7. 02152 
0. 000879 
0. 997595 
a 002243 
0. 000015 
0. 108204 
-3. 63246 
60220. 6 
643. 648 
-110. 771 
-0. 466533 
-1. 00428 
19. 6630 

19. 4297 

19. 1964 

-0.023740 
—a 44009 

-I. 47116 

^ - 

0. 030 
19. 1964 
0. 608174 
■ 0. 369876 
0. 224949 
0. 136808 ' 
81602. 9 
4. 68030 
-14. 8996 
a 001775 
2. 01448 
0. 003821- 
0. 000030 
0. 216409 
-7. 88831 
80415. 6 
555. 901 
-144. 668 
-0. 638573 
-1. 64285 
19. 0245 

18 7052 

18 3859 

-0.051552 
-4. 49248 

-1.93306. 

0. 035 
18 3859 
0. 701700 
0. 492383 
0. 346605 
0. 242441 
101248 0 
8 11688 
-26. 8749 
0. 003201 
3. 63288 
0. 005977 
0. 000054 
0. 389525 
-14. 7356 
99028 0 
572. 946 
-173. 112 
-0. 794426 
-2. 43728 
18 2300 

17. 8328 

17. 4355 

-0. 096297 
-6. 37616 

-2. 33323 
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APPENDIX B 

RELATIONS BETWEEN TWO-MASS SYSTEM AND REPRESENTED STRUCTURAL MODE 


Tile sum of the masses m£ and ms sho^ra in figure 1 must 
equal the gross mass of the represented airplane in order to 
obtain the proper nodal acceleration. For the hypothetical 
limit condition in which the wing mass is concentrated at a 
single point in each semispan, is the actual hull mass and 
ms is the actual wing mass. In order to take into account 
the more complex nature of the structural action for a par- 
ticular mode, the dete rmin ation of the ratio of these masses 
is necessary so that the vibrational energy of the simplified 
and represented sj-stems are equal for the same vibrational 
amplitudes of mt and the actual hull or point of force 
application. 

On the basis of the theory of vibrations, for which equa- 
tions are included in reference 3. the vibration energy E of 
the two-mass sj’stem is given by the equation 

ms<ps^ (BI) 

where 


Equality of the vibrational amplitude of the lower mass of 
the sLmplifieii system to the vibrational amplitude of the 
hull or fuselage of the fliung boat or airplane relative to the 
nodal point of the represented mode is expressed by the 
equation 

<PL = 'Ph CBS) 

where 

deflection of hull or fuselage of flying boat or airplane 
relative to nodal point of represented mode 

The requirement of equal energy of the simplified and 
represented sj'stems for the condition expressed by equation 
(B5) gives combination of equatipns (B3), (B4), and (B5) 
to obtain the following equation for the mass ratio of the 
two-mass system: 

m,s m ip\ . CB61 


vibrational amplitude of relative to nodal point of 
sj’stem 

<Ps vibrational amplitude of ms relative to nodal point of 
system 

w natural frequency 

Since for the two-mass system the node is at the center of 
gravity, 

mt<pL=ms9s (B2) 


If equation (Bl) and equation (B2) are combined, and since 
the total mass m is equal to the sum of the masses ms and 
fflt, the following equation can be obtained for the vibrational 
eneigj' E of the two-mass system; 


jS'_l mm£ j 
w-~2 ms 


(B3) 


The vibrational energy of the represented mode is a function 
of the spanwise mass distribution and mode shape. On the 
basis that the semispan of the airplane is divided into j sec- 
tions or stations, the vibrational eneigy of the mode can be 
written as follows: 


Computation of the natural frequency of wing modes has 
received a great deal of attention in connection with study 
of ^ving flutter and need not be treated herein. Incidental 
to calculation of the natural frequency, a mode shape is 
attained which, together with knowledge of the mass distri- 
bution, permits use of equation (B6). In cases in which the 
wing has been constructed, the mode shape and natural 
frequency may be determined experimentally. Equation 
(A4) in appendix A of the present report permits computa- 
tion of the spring constant which for a given mass ratio of 
the simplified system gives the required natural frequency'. 

After the accelerations of the two masses of the simplified 
system have been computed, equation (A3) of appendix A 
fixes the magnitude of the nodal acceleration. The diderence 
between the nodal acceleration and the huU, fuselage, or 
float acceleration can be taken aS a measure of the osdUatory' 
acceleration. On the basis that the structure deflects in the 
mode used in determining the equivalent two-mass system, 
the acceleration at any point is given by' the equation 

(B7) 

<Pk 


E 1 

—,=^ {rni(pi*+m2(p-*+rni<p3^+ . . . +mjcpj*) (B4) 

03“ Z 

where 

mj mass of jth spanwise section; value for semispan doubled 
to represent the entire span 

deflection of mass at yth spanw'ise section relative to 
nodal point (hf+xa.f) 

■ kj deflection of elastic axis at jth spanwise deflection rela- 
tive to nodal point 

X chordwise distance from elastic axis to effective mass 
center 

aj torsional deflection at jth spanwise station 


The foregoing equations, with consideration of rotatory’ 
inertia and energy, may be applied to tlie case of a tip float 
attaclied to a fle.xible wing if the stiffness and mass distribu- 
tion are known and if a manner of structural deflection is 
assumed. 
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